Ion exchange phase transitions in "doped" water— filled channels 
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Ion transport through narrow water-fiUed channels is impeded by a high electrostatic barrier. The 
latter originates from the large ratio of the dielectric constants of the water and a surrounding media. 
We show that "doping", i.e. immobile charges attached to the walls of the channel, substantially 
reduces the barrier. This explains why most of the biological ion channels are "doped". We show 
that at rather generic conditions the channels may undergo ion exchange phase transitions (typically 
of the first order). Upon such a transition a finite latent concentration of ions may either enter 
or leave the channel, or be exchanged between the ions of different valences. We discuss possible 
implications of these transitions for the Ca-vs.-Na selectivity of biological Ca channels. We also show 
that transport of divalent Ca ions is assisted by their fractionalization into two separate excitations. 
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I. INTRODUCTION 

Protein ion channels functioning in biological lipid 
membranes is a major frontier of biophysics An ion 
channel can be inserted in an artificial membrane in vitro 
and studied with physical methods. For example, one can 
measure a current-voltage response of a single water- 
filled channel connecting two water reservoirs (Fig. ^ 
as function of concentration of various salts in the bulk. 
It is well known24i^ that a neutral channel creates a 
high electrostatic self-energy barrier for ion transport. 
The reason for this phenomena lies in the high ratio of 
the dielectric constants of water, ki ~ 80, and the sur- 
rounding lipids, K2 — 2. For narrow channels the barrier 
significantly exceeds ksT and thus constitutes a serious 
impediment for ion transport across the membrane. It 
is a fascinating problem to understand mechanisms "em- 
ployed" by nature to overcome the barrier. 

At a large concentration of salt in the surrounding wa- 
ter the barrier can be suppressed by screening. How- 
ever, for biological salt concentrations and narrow chan- 
nels the screening is too weak for tha^. As a result, 
at the ambient salt concentrations even the screened 
barrier is usually well above fcsT. What seems to be 
the "mechanism of choice" in narrow protein channels is 
doping" . Namely, there is a number of amino- acids con- 
taining charged radicals. Amino-acids with, say, negative 
charge are placed along the inner walls of the channels. 
The static charged radicals are neutralized by the mobile 
cations coming from the water solution. This provides 
a necessary high concentration of mobile ions within the 
channel to suppress the barrier—. 

Similar physics is at work in some artificial devices. 
For example, water filled nanopores are studied in silicon 
or silicon oxide filma^. Dielectric constant of silicon oxide 
is close to 4 <C 80, so a very narrow and long artificial 
channels may have large self-energy barrier. Ion trans- 
port through such a channel can be facilitated by natu- 
rally appearing or intentionally introduced wall charges, 
" dopants" . Their concentration may be tuned by pH of 
the solution. 

The aim of this paper is to point out that the dop- 



ing may lead to another remarkable phenomena: the ion 
exchange phase transitions. An example of such a tran- 
sition is provided by a negatively doped channel in the 
solution of monovalent and divalent cations. At small 
concentration of divalent cations, every dopant is neu- 
tralized by a single monovalent cation. If the concen- 
tration of divalent salt increases above a certain critical 
concentration the monovalent cations leave the channel, 
while divalent ones enter to preserve the charge neutral- 
ity. Since neutralization with the divalent ions requires 
only half as many ions, it may be carried out with lesser 
entropy loss than the monovalent neutralization. The 
specifics of the Id geometry is that this competition leads 
to the first order phase transition rather than a crossover 
(as is the case for neutralizing 2d charged surfaces in the 
solution). We show that the doped channels exhibit rich 
phase diagrams in the space of salt and dopant concen- 
trations. 
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FIG. 1: Electric field of a cation in a cylindrical channel with 
the large dielectric constant ki K2. L is the channel length, 
a is its radius. The self-energy barrier is shown as a function 
of the X coordinate. 

Let us remind the origin of the electrostatic self-energy 
barrier. Consider a single cation placed in the middle of 
the channel with the length L and the radius a, Fig. ^ 
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If the ratio of the dielectric constants is large ki/k2 ^ 1, 
the electric displacement D is confined within the chan- 
nel. As a result, the electric field lines of a charge are 
forced to propagate within the channel until its mouth. 
According to the Gauss theorem the electric field at a 
distance x > a from the cation is uniform and is given by 
Eq = 2e/(Kia^). The energy of such a field in the volume 
of the channel is: 

5 2 ^ ~A ' (-^^ 

ovr ZKifl^ 4 

where the zero argument is added to indicate that there 
are no other charges in the channel. The bare barrier, 
C/i(0), is proportional to L and (for a narrow channel) 
can be much larger than fc^T, making the channel resis- 
tance exponentially large. 
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FIG. 2: A cation (in thin circle) bound to a negative wall 
charge (in thick circle). When the cation moves away from 
the host the energy grows linearly with the separation |a;|. 

If a dopant with the unit negative charge is attached to 
the inner wall of the channel it attracts a mobile cation 
from the salt solution (Fig. [21). There is a confining inter- 
action potential $(a;) — Eq\x\ between them. Condition 
e^{xT) — ksT defines the characteristic thermal length 
of such classical "atom", xt — ksT/eEQ — 121b, 
where Ib = I {i^ikBT) is the Bjerrum length (for water 
at the room temperature Ib = 0.7 nm). This "atom" is 
similar to an acceptor in a semiconductor (the classical 
length Xt plays the role of the effective acceptor Bohr 
radius). It is convenient to measure the one-dimensional 
concentrations of both mobile salt and dopants in units of 
\/xt- In such units the small concentration corresponds 
to non-overlapping neutral pairs ("atoms"), while the 
large one corresponds to the dense plasma of mobile and 
immobile charges. These phases are similar to lightly and 
heavily doped p-type semiconductors, respectively. 

It is important to notice that in both limits all the 
charges interact with each other through the Id Coulomb 
potential 

<^{xi- Xj) ^ aiajE,^\xi- Xj\, (2) 

where Xi and ct^ = ±1 are coordinates and charges of 
both dissociated ions and dopants. Another way to for- 
mulate the same statement is to notice that the electric 
field experiences the jump of 2i?oO'i at the location of the 
charge ai. Because all the charges inside the channel are 
integers in unit of e, the electric field is conserved modulo 
2i?o. It is thus convenient to define the order parameter 
q = frac[£'(a;)/2_Eo], which is the same at every point 
along the channel. The physical meaning of q G [0, 1] 



is the image charge induced in the bulk solution to ter- 
minate the electric field lines leaving the channel. One 
may notice that the adiabatic transport of a unit charge 
across the channel is always associated with q spanning 
the interval [0, 1]. Indeed, a charge at a distance x from 
one end of the channel produces the fields 2Eqx/L and 
2Eq{x/L — 1) to the right and left of x, correspondingly. 
Therefore q — x/ L continuously spans the interval [0, 1] 
as the charge moves from a; = to x = L. 

To calculate the transport barrier (as well as the ther- 
modynamics) of the channel one needs to know the free 
energy Fq of the channel as a function of the order pa- 
rameter. The equilibrium (referred below as the ground 
state) free energy corresponds to the minimum of this 
function i^min- Transport of charge and thus varying q 
within [0, 1] interval is associated with passing through 
the maximum of the Fq function i^max- Throughout this 
paper we shall refer to such a maxima as the saddle point 
state. The transport barrier is given by the difference 
between the saddle point and the ground state free en- 
ergies: Ul = ^max — -Fmin. The equilibrium concentra- 
tions of ions inside the channel are given by the deriva- 
tives of i^min with respect to the corresponding chemi- 
cal potentials related to concentrations in the bulk solu- 
tion. We show below that the calculation of the partition 
function of the channel may be mapped on a fictitious 
quantum mechanical problem with the periodic poten- 
tial. The function Fq plays the role of the lowest Bloch 
band, where q is mapped onto the quasi-momentum. As 
a result, the entire information of the thermodynamical 
as well as transport properties of the channel may be ob- 
tained from the analytical or numerical diagonalization 
of the proper "quantum" operator. 

For an infinitely long channel with the long range in- 
teraction potential Eq. (|2Jl we arrive at true phase tran- 
sitions, in spite of the one-dimensional nature of the 
problem. However, at finite ratio of dielectric constants 
electric field lines exit from the channel at the distance 
^ ~ a{Ki/K2)^^^ ~ 6.8a. As a result, the potential 
Eq. ^ is truncated and phase transitions are smeared 
by fluctuations even in the infinite channel. In practice 
all the channels have finite length which leads to an ad- 
ditional smearing. We shall discuss sharpness of such 
smeared transitions below. 

The outline of this paper is as follows: In section UTI 
we briefly review results for the simplest model of peri- 
odically placed negative charges in the monovalent salt 
solution (Fig. O, published earlier in the short communi- 
cation^. In this example the barrier is gradually reduced 
with the increasing doping. Sections IIIII - IVl are devoted 
to several modifications of the model which, contrary to 
expectations raised by the results of section ^ lead to 
ion exchange phase transitions. In section IIIII we con- 
sider a channel with the alternating positive and nega- 
tive dopants in monovalent salt solution and study the 
phase transition at which mobile ions leave the channel. 
In section IIVI and we return to equidistant negative 
dopants, but consider the role of divalent cations in the 
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FIG. 3: The ground state and the transport saddle point of 
the channel with negative monovalent dopants. Only the right 
part of the channel is shown, (a) The ground state {q — 0): 
all dopants (thick circles) bound mobile cations (thin circles), 
(b) The transport saddle point {q — 1/2): cations are free in 
sections between two adjacent dopants. 



bulk solution. In particular, in section IIVI we assume 
that all cations in the bulk solution are divalent and show 
that this does not lead to four times increase of the self- 
energy barrier. The reason is that the divalent ions are 
effectively fractionalized in two monovalent excitations. 
In section we deal with a mixture of monovalent and 
divalent cations and study their exchange phase transi- 
tion. We discuss possible implications of this transition 
for understanding the Ca-vs.-Na selectivity of biological 
Ca channels. For all the examples we present transport 
barrier, latent ion concentration and phase diagram along 
with the simple estimates, explaining the observed phe- 
nomenology. The details of the mapping onto the effec- 
tive quantum mechanics as well as of the ensuing numer- 
ical scheme are delegated to sections IVll and IVIll Results 
of sections IIIII - 13 are valid only for very long channels 
and true Id Coulomb potential. In section rVlIII we dis- 
cuss the effects of finite channel length and electric field 
leakage from the channel on smearing of the phase tran- 
sitions. In section Hxl we consider boundary effects at the 
channel ends leading to an additional contact (Donnan) 
potential. We conclude in section H by brief discussion 
of possible nano-engineering applications of the presented 
models. 



II. NEGATIVELY DOPED CHANNEL IN A 
MONOVALENT SOLUTION 

As the simplest example of a doped channel we con- 
sider a channel with negative unit-charge dopants pe- 
riodically attached to the inner walls at distance xt/"/ 
from each other (Fig. Here 7 is the dimensionless 
one-dimensional concentration of dopants. At both ends 
(mouths) the channel is in equilibrium with a monovalent 
salt solution with the bulk concentration c. It is conve- 
nient to introduce the dimensionless monovalent salt con- 
centration as ai = cko^xt- We shall restrict ourselves 
to the small salt concentration (or narrow channels) such 
that q;i ^ 1. In this case the transport barrier of the 
undoped (7 = 0) channel is given by Eq. (save for the 



small screening reduction which scales as 1 — 4ai, ja])- 

The calculations, described in details in sections IVII 
and IVIII lead to the barrier plotted in Fig. 0] as a func- 
tion of the dopant concentration 7. The barrier decreases 
sharply as 7 increases. For example, a very modest con- 
centration of dopants 7 — 0.2 is enough to suppress the 
barrier more than five times (typically bringing it below 
ksT). There are no phase transitions in this system in 
the entire phase space of concentrations 7 and ai. 
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FIG. 4: The function I7l(7)/I7l(0) for Qi = 10"^ Its 7 < 1 
asymptotic, Eq. 101 , is shown by the dotted line. 

For the small dopant concentration, 7^1, the result 
may be understood with a simple reasoning. The ground 
state of the channel corresponds to all negative dopants 
being locally neutralized by the mobile cations from the 
solution (Fig. Ola). As a result, there is no electric field 
within the channel and thus in the ground state q = 0. 
The corresponding free energy has only entropic compo- 
nent (the absence of the electric fields implies no energy 
cost) which is given by: 

Fo = -7LfcsTln(2ai) = Ul{0) ■ 47ln[l/(2ai)] . (3) 

Indeed, bringing one cation from the bulk solution into 
the channel to compensate a dopant charge leads to the 
entropy reduction Sq = fc^ ln(7ra^2a;Tc) = fcBln(2Q;i). 
Here 7ra^2xT is the allowed volume of the cation's ther- 
mal motion within the channel, while is the accessible 
volume in the bulk. 

The maximum of the free energy is associated with 
the state with q = 1/2, see Fig. ^h. It can be viewed 
as a result of putting a vacancy in the middle of the 
channel. The latter creates the electric field zLEq which 
orients the dipole moments of all the "atoms" . In other 
words, it orders all the charges in an alternating sequence 
of positive and negative ones. This unbinds the cations 
from the dopants and makes them free to move between 
neighboring dopants (Fig.|31b). Indeed, upon such rear- 
rangement the electric field is still E = ±Eo everywhere 
in the channel, according to the Gauss theorem. There- 
fore, the energy of g = 1/2 state is still given by Eq. 
However, its entropy is dramatically increased with re- 
spect to g = state due to the unbinding of the cations: 
the available volume is now tto^xt/'J and the resulting 
entropy per cation is S'1/2 — kghi^ai/j). The corre- 
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spending free energy of the saddle point state is: 



Ut 



-Fi/2 = t^L(0)[l-47ln(ai/7)] 



(4) 



Recalling that the transport barrier is given by the dif- 
ference between the saddle point and the ground state 
free energies, one obtains: 



C/L(7) = t/L(0)[l-47ln(l/27)] 



(5) 



This expression is plotted in Fig.0]by the dotted line. It 
provides a perfect fit for the transport barrier at small 
dopant concentration. Equation lO is applicable for 
ai < 7 <C 1. In the opposite limit 7 < ai <C 1 more 
free ions may enter the channel in the saddle point state. 
As a result, the calculation of 51/2 should be slightly 
modified^, leading to: 



C/l(7) = Ul{Q) 



1 (X\ 

1 — 47 In ( sinh — 
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(6) 



This result, exhibiting non-singular behavior in the small 
concentration limit, is valid for an arbitrary relation be- 
tween ai and 7 (both being small enough). 
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FIG. 6: The transport barrier in units of Ul{Q) for the "com- 
pensated" channel of Fig. |3 for cti = 0.01. The dotted line 
shows Eq. 1)10^ . 



at q — 1/2, i.e. the electric field ±£'9. However, once 
the dopant concentration 7 increases the second mini- 
mum develops at q — 1/2, which eventually overcomes 
the q — Q minimum at 7 = 7^. In the limit of large 7 the 
ground state corresponds to g = 1/2 (electric field ±£^0)7 
while the saddle point state is at q = (no electric field). 
It is clear from Fig.[7|that the transition between the two 
limits is of the first order. 



III. CHANNEL WITH ALTERNATING 
POSITIVE AND NEGATIVE DOPANTS 

As a first example exhibiting the ion-exchange (actu- 
ally ion-release) phase-transition we consider a model 
of a "compensated" channel (the word compensation is 
used here by analogy with semiconductors, where accep- 
tors can be compensated by donors) . This is the channel 
with positive and negative unit-charge dopants alternat- 
ing in one-dimensional NaCl type lattice with the lattice 
constant 2x^/7 (Fig. 13 . 
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FIG. 5: A channel with the alternating dopants (thick circles). 
The mobile counterions are shown as thin circles. 

The channel is filled with the solution of the mono- 
valent salt with the bulk dimensionless concentration 
ai = tto^xtc. The transport barrier calculated for 
Qfi = 0.01 as a function of the dopant concentration 
7 is depicted in Fig. |B1 One observes a characteristic 
sharp dip in a vicinity of a certain dopant concentration 
7c ~ 0.06. To clarify a reason for such an unexpected 
behavior (cf. Fig. ^ we plot the free energy as a func- 
tion of the order parameter q for a few values of 7 close 
to 7c, see Fig.Q Notice that for small 7 the minimum of 
the free energy is at q = 0, corresponding to the absence 
of the electric field inside the channel. The maximum is 
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FIG. 7: Free energy in units of Ul{0) as a function of q for 
7 = 0.02, 0.05, 0.07 and 0.09 (from top to bottom) at ai = 
0.01. The lower three graphs are vertically offset for clarity. 

To understand the nature of this transition, consider 
two candidates for the ground state. The q — state, 
referred below as 0, is depicted on Fig.|Sla. In this state 
every dopant tightly binds a counterion from the solution. 
Such a state does not involve an energy cost and has 
the negative entropy Sq — fcBln(2ai) per dopant. As a 
result, the corresponding free energy is (cf. Eq. Q): 

Fo-C/L(0)-47ln[l/(2ai)]. (7) 

An alternative ground state is that of the channel free 
from any dissociated ions. Fig. |S1 b. There is an elec- 
tric field zLEq alternating between the dopants. This is 
g = 1/2 state, or simply 1/2 state. Since no mobile ions 
enter the channel, there is no entropy lost in comparison 
with the bulk. There is, however, energy cost for having 
electric field ±Eq across the entire channel. As a result, 
the free energy of the 1/2 state is (cf. Eq. (Q): 

Fi/2 = ULiO). (8) 
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Comparing Eqs. |7J and ((SJ, one expects that the critical 
dopant concentration is given by 

7c = -[41n(2al)]-^ (9) 

For 7 < 7c the state is expected to be the ground state, 
while for 7 > 7c the state 1/2 is preferable. 
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FIG. 8: States 0, 1/2 and 00. The corresponding free energies 
are given by Eqs. 0, ||SJ and lITTll . 



Crossing the phase boundary in either direction is as- 
sociated with an abrupt change in the one-dimensional 
density of the salt ions within the channel. The latter 
may be evaluated as the derivative of the ground state 
free energy with respect to the chemical potential of the 
salt: rijQj^ = —j^^^dFj^in/dai. In Fig. ^| we plot con- 
centration of ions within the channel, tliq^j, in units of 
dopant concentration j/xt as a function of the bulk salt 
normality ai. One clearly observes the latent concen- 
tration associated with the first-order transition. As the 
bulk concentration increases past the critical one, the mo- 
bile ions abruptly enter the channel. One can monitor the 
latent concentration An along the phase transition line 
of Fig. In Fig. ^] we plot the latent concentration 
along the phase boundary as function of the critical ai. 
As expected, in the dilute limit the latent concentration 
coincides with the concentration of dopants (i.e. every 
dopant brings one mobile ion). On the other hand, in 
the dense limit, the latent concentration is exponentially 
small (but always finite). 



In Fig. Owe plot the phase diagram on the (7, ai) 
plane. The phase boundary between and 1/2 states is 
determined from the condition of having two degenerate 
minima of the free energy Fg a.t q = and q — 1/2. In 
the small concentration limit (see the inset in Fig.l^J the 
phase boundary is indeed perfectly fitted by Eq. Q . For 
larger concentrations it crosses over to 7c oc ^/al. This 
can be understood as a result of a competition between 
dopant separation xt/j and the Debye screening length 
r/3 cx xt/^/oi- For 7 < ^/ai we have ru < xt/j, and 
thus each dopant is screened locally by a cloud of mobile 
ions. As a result, the neutral state is likely to be the 
ground state. In the opposite case 7 > ^/al the Debye 
length is larger than the separation between oppositely 
charged dopants. This is the limit of a week screening 
when most dopants may not have counterions to screen 
them. Thus the 1/2 state has a chance to have a lower 
free energy. 
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FIG. 9: The phase diagram of the channel with alternating 
doping (dotted lines). The phase boundary between the q — 
and q = 1/2 phases can be fitted as ai ~ 67^ for ai > 1. 
Inset: the small concentration part of the phase boundary 
fitted by Eq. @ (the full line). 
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FIG. 10: The concentration of cations inside the channel in 
units of y/xT for 7 = 0.1. The discontinuous change occurs 
at the phase transition point. 
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FIG. 11: The latent concentration of cations in units oi y/xr 
along the phase boundary line. 

Let us return to the calculation of the transport bar- 
rier. Fig. El To this end one needs to understand the na- 
ture of the saddle point states (in addition to that of the 
ground states, discussed above). Deep in the phase where 
is the ground state, the role of the saddle point state is 
played by the 1/2 state. Correspondingly the transport 
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barrier is approximately given by the difference between 
Eq. ||SJ) and Eq. {Tj). One may improve this estimate by 
taking into account that in the 1/2 state the free ions 
may enter the channel (in even numbers to preserve the 
total charge neutrality). This leads to the entropy of the 
1/2 state given by Sy^ = fcs lnEr=o(«i/7)'V(2fc)!] per 



dopant. As a result, 
in the state: 

UL{ai,7) = ULiO) 



one finds for the transport barrier 
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(10) 



This estimate is plotted in Fig. El by the dotted line. 



IV. NEGATIVELY DOPED CHANNEL WITH 
DIVALENT CATIONS 

In the above sections all the mobile ions as well as 
dopants were monovalent. In this section we study the 
effect of cations being divalent (e.g. Ca^^ , or Ba^+), 
while all negative charges (both anions and dopants) are 
monovalent (Fig. I13|l . For example, one can imagine a 
channel with negative wall charges in CaCl2 solution. We 
denote the dimensionless concentration of the divalent 
cations as a2- The concentration of monovalent anions 
is simply a-i — 1oi2- 
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FIG. 12: Free energy diagram for a\ = 0.01 as function of 
dopant concentration 7. The solid lines are numerical results 
for the ground and saddle point states. Eqs. (O, l|SJ and lllL 
describing states 0, 1/2 and 00, correspondingly, are shown 
by dashed, dotted and dash-dotted lines. 

In the 1/2 state, Fig.|Slb, the channel is almost empty 
in the ground state configuration. The saddle point may 
be achieved by putting a single cation in the middle of 
the channel. This will rearrange the pattern of the in- 
ternal electric field as depicted in Fig. [S] c. This state 
corresponds to q = and is denoted as 00 (to distinguish 
it from the state 0, Fig.|Hla). Its free energy (coinciding 
with the energy) is given by: 



^^oo = 2C/i(0) 



(11) 



In Fig. El we plot the free energies of the three states 0, 
1/2 and 00 (cf. Fig.|Hland Eqs. ©, © and IJlJ) as func- 
tions of the dopant concentration 7. On the same graph 
we also plot the calculated ground state free energy F^in 
along with the saddle point free energy F^ax- It is clear 
that the ground state undergoes the first order transition 
between and 1/2 states at 7 = 7c. On the other hand, 
the saddle point state experiences two smooth crossovers: 
first between 1/2 and and second between and 00 
states. The difference between the saddle point and the 
ground state free energies is the transport barrier, which 
exhibits exactly the type of behavior observed in Fig. 
Curiously, at large concentration of dopants the barrier 
approaches exactlythe same value as for the undoped 
channel, [/^(ai), p. This could be expected, since ex- 
tremely closely packed alternative dopants compensate 
each other, effectively restoring the undoped situation. 
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FIG. 13: A channel with periodic dopants (thick circles) and 
divalent cations. Mobile ions (Ca2+ and CP) are in thin 
circles. 

The transport barrier as function of the dopant concen- 
tration 7 for a2 = 5 • 10~^ is shown in Fig. 1141 Similarly 



0.5 - 




Yc 0.05 0.1 



Y 



FIG. 14: The transport barrier [/L(a2,7) in units of [/^(O) 
for 02 = 5 • 10~^ (sohd line). The dotted line is 1-^1/2 — fb|, 
and dashed line is |_Foo —^^1/2 1 calculated using Eqs. I12L I13II 
and 11511 . 

to the case of the alternating doping (cf. Fig. (BJ the 
barrier experiences a sudden dip at some critical dopant 
concentration 7c ~ 10~^. To understand this behavior 
we looked at the free energy Fq as function of the or- 
der parameter q for several values of 7 in the vicinity of 
7c. The result is qualitatively similar to that depicted in 
Fig. El Thus, it is again the first order transition between 
two competing states that is responsible for the behavior 
observed in Fig. ITll 

The two candidates for the ground state, denoted, in 
accordance with the fractional part of the internal electric 
field g, as and 1/2, are depicted in Fig. E|a,b. In the 
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FIG. 15: The states 0, 1/2, and 00 for the channel with di- 
valent cations. The corresponding free energies are given by 
Eqs. I12II . 1131 and 1151 1. Dopants (thick circles) and anions 
(thin circles) are monovalent, while cations (shown by 2+ in 
thin circles) are divalent. Electric fields are shown schemati- 
cally, one line per Eq. 



evolves smoothly from 1/2 to and eventually to 00 (up- 
per full line). The difference between the two gives the 
transport barrier depicted in Fig. where \Fi^2 ^ Po\ 
and ji^oo ~ ^1/2! are also shown for comparison. 
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state every dopant is screened locally by an anion and 
a doubly charged cation (Fig. E| a) . The corresponding 
free energy (consisting solely from the entropic part) is 
given b}*^: 



Fo = C/.(0).4,ln(^ 



(12) 



The other state, 1/2, has every second dopant over- 
screened by a divalent cation (Fig. E]b) . The free energy 
of this state is: 



Fi/2 = C/l(0) 



l + |^^ln(l/a2) 



(13) 



Comparing Eqs. (|12|l and (|13|l . one finds for the critical 
dopant concentration: 



7c = - [61n(Q2) +41n6]" 



(14) 



For a2 = 5 • 10~^ it leads to 7c ~ 0.012 in a good agree- 
ment with Fig. 

To explain the transport barrier observed for small 7 
(Fig. [T^ one needs to know the saddle point state. Such 
a state is depicted in Fig.^jc and is denoted as 00. It 
has one divalent cation trapped between every other pair 
of dopants. It is easy to see that for such arrangement 
the cations are free to move within the "cage" defined by 
the two neighboring dopants. This renders a rather large 
entropy of the 00 state. Its free energy is given by 



f/L(0) 



2 + Yln(7/a2) 



(15) 



In Fig. ^] we plot the free energies of the states 0, 1/2 
and 00, given by Eqs. H12|l . (|13l) and (|15|l correspond- 
ingly, as functions of 7. On the same graph we plot 
calculated ground state free energy Fmin along with the 
saddle point free energy i^max (full lines). One observes 
that the ground state indeed undergoes the first order 
transition between the states and 1/2 upon increasing 
7 (lower full line). On the other hand, the saddle point 



FIG. 16; Free energy diagram for Q2 = 5 • 10"*^ as a function 
of dopant concentration 7. The full lines are numerical results 
for the ground and saddle point states respectively. Eqs. I13I I. 
(I15II and 11211 describing states 1/2, 00 and 0, are shown by 
dashed, dash-dotted and solid thin line correspondingly. State 
is the ground state for 7 < 7c, state 1/2 is the ground state 
for 7 > 7c, and the state 00 is the saddle point for 7 > 7c. 
The ground state undergoes the first order phase transition at 
7c, while the saddle point state evolves in a continuous way. 

It is worth noticing that on the both sides of the tran- 
sition the transport barrier is close to Ul{0), character- 
istic for the transfer of the unit charge e. One could 
expect that for charges 2e the self-energy barrier should 
be rather 4f7i(0). This apparent reduction of the charge 
seems natural for the very small 7. Indeed, the corre- 
sponding ground state (Fig. I15| a) contains complex ions 
(CaCl)^+ with the charge e. On the other hand, for 
7 > 7c the channel is free from Cl~ ions and the cur- 
rent is provided by Ca^"*" ions only. Thus, the observed 
barrier of < Ul (0) may be explained by fractionalization 
of charges 2e into two charges e. To make such a frac- 
tionalization more transparent one can redraw the saddle 
point state 00 in a somewhat different way, namely cre- 
ating a soliton (domain wall, or defect) in the state 1/2. 

Fig. El shows the channel with a soliton in the mid- 
dle. One can see that in this version of the 00 state the 



FIG. 17: The saddle point state 00 represented as the unit 
charge defect (soliton) within the ground state 1/2. Ca'^'*' ion 
added to the channel in the state 1/2 fractionalizes in two 
such solitons, each having an excess energy < Ul{0). 

fields 2Eq and alternate similarly to Fig. El c. We can 

also see that the soliton in the middle has the charge 
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e. Fractionalization of a single Ca^+ ion in two charge- 
e defects means that a Ca ion traverses the channel by 
means of two solitons moving consecutively across the 
channel. The self-energy of each soliton (and therefore 
the transport barrier) does not exceed Ul{0)- 
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V. NEGATIVELY DOPED CHANNEL IN 
SOLUTION WITH MONOVALENT AND 
DIVALENT CATIONS 

We turn now to the study of the channel in a solu- 
tion with the mixture of monovalent cations with the di- 
mensionless concentration ai and divalent cations with 
the concentration a2- Neutrality of the solution is 
maintained by monovalent anions with the concentration 
a_i = ai + 2a2- The channel is assumed to be doped 
with the unit charge negative dopants, attached periodi- 
cally with the concentration 7. 

In Fig. 1201 we plot the barrier as a function of the diva- 
lent cation concentration a2 for 7 — 0.1 and ai — 0.001. 
The overall decrease of the barrier with the growing ion 
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FIG. 18: ln[[/L(7)/C/z,(0)] for Q2 = 5* 10"^ Each dip signals 
the presence of a phase transition. 
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FIG. 19: The phase diagram of the channel with the divalent 
cations on the (7, 02) plane. The dotted lines are the bound- 
aries between phases. The ground states in each region are 
labelled. Inset: magnification of the first phase boundary line 
at small 7. The solid line is Eq. H14|l . 

So far the physics of the first order transition in the 
divalent system was rather similar to that in the channel 
with the alternating doping and a monovalent salt. There 
are, however, important distinctions, taking place at 
larger dopant concentration. The logarithm of the trans- 
port barrier in the wide range of 7 is plotted in Fig. 1181 
One notices a series of additional dips at jci ~ 1-7, 
7c2 ~ 5.3, 7c3 ~ 10.8, etc. (the first order transition at 
7c ~ 0.01, discussed above, is not visible at this scale). 
These dips are indications of the sequence of reentrant 
phase transitions, taking place at larger 7. The calcu- 
lated phase diagram is plotted in Fig. E| The leftmost 
phase boundary line corresponds to the first order phase 
transition between and 1/2 states. Its low concentra- 
tion part along with the fit with Eq. H14|) is magnified 
in the inset. The other lines are transitions spotted in 
Fig. (|18|l . We discuss them in Appendix IXI 
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The transport barrier as a function of In 0:2 for 7 
0.1 and ai = 0.001. 



FIG. 20: 



concentration is interrupted by the two sharp dips at a2c 
and a'2c- By plotting the Fq function for several a2 in the 
vicinity of a2c and a'2c, one observes that they correspond 
to the two consecutive first order transitions. As a2 in- 
creases the system goes from q ~ phase into q = 1/2 
phase and eventually back into q = phase. 
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FIG. 21: Concentration of Ca^+ ions in the channel in units of 
y/xT as function of In 02 for 7 = 0.1 and Qi = 0.001. There 
are two discontinuous changes in Ca concentration at a2c and 
Qf2c- Each of these two latent changes is close to the half of 
the number of dopants. 
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The concentration of divalent cations within the chan- 
nel is given by 



nca = 



0L2 



dF,r 



knTL dai 



(16) 



It is plotted as a function of the bulk concentration in 
Fig. 1^ One observes that at the first transition the 
number of divalent cations entering the channel is close 
to the half of the number of dopants. When the second 
transition is completed the number of divalent cations is 
approximately the same as the number of dopants. This 
provides a clue on the nature of the corresponding ground 
states. For a2 < a2c there are almost no divalent cations 
in the channel. Therefore, both the ground state and the 
saddle point state are the same as in sec. m shown in 
Fig. 13 The ground state free energy and the transport 
barrier are given by Eqs. Q and (|SJl correspondingly. At 
Q!2 = a2c the first order ion-exchange phase transition 
takes place, where every two monovalent cations are get- 
ting substituted by a single divalent one. The system's 
behavior at larger 012 is qualitatively similar to that de- 
scribed in the previous section. For a2c < a2 < a'2c 
the ground state is the 1/2 state, pictured in Fig.lTslb. 
The corresponding free energy is given by Eq. (|13|l . The 
second phase transition at a2c is similar to that taking 
place on the leftmost phase boundary of Fig. ^| (which 
is crossed now in the vertical direction). For a2 > 
the ground state is the state (Fig. E| a) with the free 
energy given by Eq. ifT^ . 



Notice that it scales as af and therefore at small con- 
centrations the transition takes place at a2c ^ cti. This 
is a manifestation of the law of mass action. The sec- 
ond critical value may be estimated from Eq. I|14|) as 
a'2^ — e"^/^^''') and is approximately independent on ai. 
Comparing the two, one finds that the transitions may 
take place only for small enough concentration of the 
monovalent ions Inai ^ —7/3 — In 2. For larger ai there 
is a smooth crossover between the two q = states. 




FIG. 23: Phase diagram in the space of cation and dopant 
concentrations (lnQ2 — 7 — Inai). Inside the tent-like shape 
the ground state is q = 1/2, outside - the ground state is 

q = 0. 
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FIG. 22: Free energies of the three competing ground states 
for 7 = 0.1 and m = 0.001. The dash-dotted line is Eq. ©, 
the solid line is Eq. 11311 . and the dashed line is Eq. 1121 1. The 
actual ground state is chosen as the lowest of them. 

The free energies of the three competing ground states 
as functions of a2 are plotted in Fig. [521 for the same 
parameters as in Fig. |2D1 They indeed intersect at the 
concentrations close to a2c and a2c- Since at each such 
intersection the symmetry of the ground state (the q- 
value) changes, one expects that the ground state changes 
via the first-order phase transition. The critical value a2c 
of the ion-exchange transition may be estimated from 
Eqs. (O and Q as: 



a2c 



(2ai)2ei/(27) 



The phase diagram in the space of cation and dopant 
concentrations (lna2 — 7 — Inai) is plotted in Fig.|2Sl By 
fixing some 7 and not too large ai and varying a2 one 
crosses the phase boundary twice. This way one observes 
two first order phase transitions: from to 1/2 and then 
from 1/2 to 00. The corresponding transport barrier is 
shown in Fig. EOl 
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FIG. 24: Schematic plot of the current through a Ca channel 
as a function of Ca^"*" concentration, cf. Ref. flB . The Na^ 
concentration is ai ~ lO"'^, the wall charge concentration is 



7 ~ 0.3, and the (smeared) transition is at In a2c 



-10. 



(17) 



The model presented in this section is a simple cartoon 
for Ca^+ selective channels. Let us consider the total cur- 



10 



rent through the channel, /, equal to the sum of sodium 
and calcium currents. Each of these partial currents in 
turn is determined by two series resistances: the channel 
resistance and the combined contact (mouth) resistances 
(the latter are inversely proportional to the concentra- 
tion of a given cation in the bulk). We assume that at 
the biological concentration of Na+ the current I = Iq 
and discuss predictions of our model of negatively doped 
channel at 7 « 0.3 regarding I/Iq with growing con- 
centration a2 of Ca^+ cations (Fig. For 02 < a2c 
the channel is populated by Na+ and Na+ current dom- 
inates. For a2 > a2c the Na"*" current is blocked because 
Na+ ions are expelled from the channel and substituted 
by the Ca^+. Indeed, a monovalent Na+ cation cannot 
unbind the divalent Ca^+ from the dopants. Therefore, 
the transport barrier for Na^ ions is basically the bare 
barrier C/l(0) S> fc^T. In these conditions the Ca^+ re- 
sistance of the channel is small, because Ca^+ concentra- 
tion in the channel is large. However, at the transition 
a2 — ot2c <C OLi the concentration of Ca^"*" ions in the 
bulk water is so small that the contact Ca^+ resistance 
is very large. Therefore, the Ca^"*" current is practically 
blocked and the total current, /, drops sharply at a2c- 
As a2 grows the Ca^"^ current increases proportional to 
a2 due to decreasing contact resistance. As a result, one 
arrives at the behavior of ///q schematically shown in 
Fig. 1241 This behavior is in a q ualita tive agreement with 
the experimental data of Refs. |lllO|| . 



VI. ANALYTICAL APPROACH 

Consider a gas consisting of N mobile monovalent 
cations and N' mobile monovalent anions along with 
non-integer boundary charges q and q' placed at a; = 
and X — L correspondingly. We also consider a single 
negative unit dopant charge attached at the point xq in- 



J 



p{x) - 



1 = jvp{x) sl^ 

= jj Vp{x)Ve{x) exp 



side the channel: Q < xq < L. The resulting charge 
density takes the form: 

N+N' 

p{x) = ajS{x — Xj) + qS{x) + q'S{x — L) — 5{x — xq) , 

(18) 

where Xj stay for coordinates of the mobile charges and 
Gj = ±1 for their charges. The interaction energy of such 
a plasma is given by: 



^=2 



dxdx' p{x)^{x — x')p{x') 



(19) 



where the Id Coulomb potential ^{x) — <i>(0) — Eo\x\ is 
the solution of the Poisson equation: = —2EqS{x). 

(The self-energy e$(0) will be eventually taken to infinity 
to enforce charge neutrality). 

We are interested in the grand-canonical partition 
function of the gas defined as 

00 



Ar,A"=0 



N+N' 




-UKksT) 



where p, is the chemical potential (the same for cations 
and anions) and Iq is a microscopic scale related to the 
bulk salt concentration as c = c'^/'^^'^/Zq. Factor -ko? /Iq 
originates from the integrations over transverse coordi- 
nates. 

To proceed with the evaluation of ZL{q,q') we intro- 
duce the resolution of unity written in the following way: 



N+N' 

E 



qe{Q)^q'e{L) + e{xQ) 



dx 9{x)p{x) 



N+N' 



Substituting this identity into Eq. 1)20(1 . one notices that the integrals over Xj decouple and can be performed in- 
dependently. The result of such integration along with the summation over iV and N' is exp{27ra^c dxcos9{x)}. 
Evaluation of the Gaussian integral over p{x) yields the exponent of d{x)^~^{x — x')0{x'). According to the Poisson 
equation the inverse potential is (f>~^(a; — x') — —{2Eq)~^5{x — a;')9|. As a result, one obtains for the partition 
function: 



Z{q,q) = 



deodOLde,, 
(27r)3 



^iq9i-\-iq 9 f ~iO 



=0 jve{x) exp l~Jdx 



XT 



2ai 

XT 



■ COS f x 
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where ai = tto^xtc. The integral over 9(x) runs over 
all functions with the boundary conditions 9(0) = Oq , 
e{L)=eL and eixo) = 0,,. 

It is easy to see that this expression represents the 
matrix element of the following T-exponent (or rather 
X-exponent) operator: 

Z{q, q') - {q\e-^ " e'^' " \q') , (21) 

where the Hamiltonian is given by H — [idgY ~ 2a;i cos 9 
and \q) is the eigenstate of the momentum operator idg. 
Since the Hamiltonian conserves the momentum up to an 
integer value, q' = q + M, one can restrict the Hilbert 
space down to the subspace with the fixed fractional part 
of the boundary charge < g < 1. In this subspace one 
can perform the gauge transformation, resulting in the 
Mathieu Hamiltonian with the "vector potential"iSiii*i^: 

^ (ide-q)^ ~2aiC0se. (22) 

It acts in the space of periodic functions: ^'(0) = ^'(0 + 
27r). Finally, taking the "democratic" sum over all integer 
parts of the boundary charge (with the fixed fractional 
part q), one obtains: 

Z{q) = Tr je"-^'^ e"** g"-^'"^ } . (23) 

In the more general situation the solution contains a 
set of ions with charges (valences) m G Z and the cor- 
responding dimensionless concentrations «,„ . The condi- 
tion of total electro- neutrality demands that: 

J2ma,n = 0- (24) 

rn 

The Mathieu Hamiltonian (|22|1 should be generalized 
asii: 

H, = {ide-q)^ -Y.a„,^^' . (25) 

m 

Despite of being non-Hermitian, this Hamiltonian still 
possesses a real band-structure^^ eg"''' . It is safe to assume 
that monovalent anions are always present in the solu- 
tion: a-i > 0. This guarantees that the band-structure 
has the unit period in q. 

Consider first an undoped channel. Its partition func- 
tion is given by Z{q) = Tr{exp(— ffgL/xr)}. In the long 

channel limit, L ^ xt, only the ground state, e^"'', of 
the Hamiltonian (|22|l or H25() contributes to the partition 
function. As a result, the free energy of the Id plasma is 
given by 

Fq = kBT€f{a,n)L/xT = 4C/L(0)ef )(a™) . (26) 

The equilibrium ground state of the plasma corresponds 
to the minimal value of the free energy. For the Mathieu 
Hamiltonian (|22(l eq"' has a single minimum aX q — Q 
(no induced charge at the boundary). It is also the case 



for the more general Hamiltonians (|25|) . Therefore the 
equilibrium state of a neutral Id Coulomb plasma does 
not have a dipole moment and possesses the reflection 
symmetry. 

The adiabatic transfer of the unit charge across the 
channel is associated with the slow change of q from q = 
to g = ±1. In this way the system must overcome 
the free energy maximum at some value of q (for the 
operators (|22|l and H25|) the maximum is at q = 1/2). 
As a result, the activation barrier for the charge transfer 
is proportional to the band-width of the lowest Bloch 
band&: 

Notice that in the ideal Id Coulomb plasma the trans- 
port barrier scales as the system size. It is also worth 
mentioning that both the equilibrium free energy = 
4[/2^(0)e||^j^(am) and the transport barrier, Eq. l(77|l . 
are smooth analytic functions of the concentrations a™. 
Since there is the unique minimum and maximum within 
the interval < q < 1, there are no phase transitions in 
the undoped channels. 

Most of biological ion channels have internal "dop- 
ing" wall charges within the channel. If integer charges 
rii,ri2, . . . ,?^Ar are fixed along the channel at the coor- 
dinates < xi < X2 < ■ . . < xn < L, the partition 
function is obtained by straightforward generalization of 
Eq. 123): 

Z{q) = Tr je-^^e^^-'e-^^ . . . e^^^e"^^ | . 

(28) 

As long as all are integer, the boundary charge q is 
a good quantum number of the operator under the trace 
sign. As a result, the partition function is again a peri- 
odic function of q with the unit period. 

For the sake of illustration we shall focus on systems 
with periodic arrangements of the wall charges. In this 
case the partition function (|28|l takes the form: Z{q) = 

Tr{lWg) }, where N is the number of dopants in the 

channel and lAq is the single-period evolution operator. 
We shall define the spectrum of this operator as: 

Uq ■^\^\e) = ■^\^\9) , (29) 

where 7 is the dimensionless concentration of dopants, 
defined as 7 = xtN/ L. The evolution operator lAq is 
non-Hermitian. Its spectrum is nevertheless real and 
symmetric function of q. The proof of this statement 
may be constructed in the same way as for the operator 
in Eq. (|25|! The free energy of a long doped system is 

given by Fq = ksTef'^Nh = AUL{0)ef'^ (cf. Eq. ^). 
The equilibrium ground state is given by the absolute 
minimum of this function. Similarly, the transport bar- 
rier is given by Eq. H27() . 

The simplest example is the periodic sequence of unit- 
charge negative dopants {xkj^i — Xk = L/N and — —1) 



12 



in the monovalent salt solution, section |nj The single- 
period evolution operator takes the form: 



-id. 



(30) 



with Hq given by Eq. H22|l . As shown in Ref.0its ground 

state ef\ai,"f) is a function qualitatively similar to the 
lowest Bloch band of the Mathieu operator in Eq. (|22|l . 
As a result, both equilibrium free energy and the trans- 
port barrier are smooth function of the salt concentration 
ai and the dopant concentration 7. 

The examples of sections llVl andlVlare described by the 
evolution operators, which have the form of Eq. (|30l) with 
the generalized Hamiltonian, Eq. 1)25(1 . In the example of 
section Hvl there are two non-zero concentrations: a_i = 
2a2, while in section^one deals with three types of ions: 
a-i = ai + 2a2- Finally, the alternating doping example 
of section IIIII is described by the evolution operator of 
the form 



^-ie -HJf +i9 -HJ-y 



(31) 



with the Hamiltonian H22|l . Such operators may have 
more complicated structure of their lowest Bloch band. 
In particular, the latter may have more than one mini- 
mum within the period of the reciprocal lattice: q € [0,1]. 
The competition between (and splitting of) the minima 
results in the first (second) order phase transitions. 

The analytic investigation of the spectrum of the evo- 
lution operators, such as Eq. ((30|l . is possible in the limits 
of small and large dopant concentrations 7. For 7^1 
it follows from Eq. (|30|l that only the lowest eigenval- 
ues of Hq are important. It is then enough to keep only 
the ground state of Hq, save for an immediate vicinity 
of q ~ 1/2, where the ground and the first excited states 
may be nearly degenerate. If also am <C 1 the spectrum 
of Hq along with the matrix elements of e^*^ may be 
calculated in the perturbation theory. Such calculations 
lead to the free energies of the ground and saddle point 
states, which are identical to those derived in sections 
HlMVl using simple energy and entropy counting. 

In the limit 7^1, a™ one may develop a variant of 
the WKB approximation in the plane of the complex 6'ji. 
It shows that the transport barrier and latent concentra- 
tion scales as exp{— cy^}, where non-universal numbers 
c are given by certain contour integrals in the complex 
plane of 6*. In the example of section |n] we succeeded in 
quantitative prediction of the coefficient c, see Ref. 0- 
In the case of the section IIVI the coefficient in front of 
cos(27r(7) (see Eq. f\Al\i ) is an oscillatory function of 7. 
This probably translates into the complex value of the 
corresponding c-constant. We did not succeed, however, 
in its analytical evaluation. 



VII. NUMERICAL CALCULATIONS 

By far the simplest way to find the spectrum of Uq is 
numerical. In the basis of the angular momentum e''^^ 



the Hamiltonian (|25|l takes the form of the matrix: 



k,k' 



{k + qfSk,k' - ^ a,nSk,k 



(32) 



In the same basis the dopant charge n creation opera- 
tor takes the matrix form: [e™^]^^,, = [<5fe,fc'+n], where 
k = . . . , —2, —1, 0, 1, 2, ... . Truncating these infinite ma- 
trices with some large cutoff, one may exponentiate the 



Hamiltonian and obtain the matrix form of 



The 



k.k' 



latter may be numerically diagonalized to find the func- 
tion e^^\arm"t)- The free energy and the transport bar- 
rier are then given by Eqs. (|26|l and H27() . 

We start from the simplest case of section ^ 
The monovalent salt with the concentration a^i = 
oi\ leads to the term — ai(5fc,fc'+i — Q!ii5fc^fe/_i in 



the matrix 



For illustration we show the 



4x4 



truncation of 

(i + g)' - 



and 



and [e ' 




matrices 






(-2 + 9) 



ai 
2 



For reasonable precision 



ai [Q + qf 
-ai 

/O 1 
10 
1 
Vo 

the truncation size of the numerical calculation has to be 
much larger. Typically we used 40 x 40 and checked that 
a further increase does not affect the results. 

To calculate the "energy" band Fq of a long channel 
for certain ai and 7 one can use the matrix form of 
Hq and e~*^, equate the largest eigenvalue of Uq in Eq. 
(EOJ to e^^^'/T^ and calculate Fq = kBTef^L/xr- Such 
a calculation gives the transport barrier C/i(Q!i,7) = 
Fmax - Fmin showu in Fig. m 

The models in sections IIVI and |3 are treated simi- 
larly, except that they have salt ion terms ~a25k,k'+2 — 
2a2i5fc,fc'_i and-ai(5fe,fc'+i - a25k,k'+2 - {oti -\- 2a2)Sk,k'-i 



respectively in Hq . For the model of section IIIII Hq 
contains salt ion term -'aiSk,k'+i ~ cti^k.k'-i- However, 
Uq is of the form of Eq. H31|) and its largest eigenvalue is 
denoted as e^^'? 



VIII. EFFECTS OF THE FINITE LENGTH AND 
THE ELECTRIC FIELD ESCAPE 

The consideration of the previous sections was cer- 
tainly an idealization that neglected several important 
phenomena. The most essential of them are: (i) the fi- 
nite length L of the channel; (ii) the escape of the electric 
field lines from the water into the media with smaller 
dielectric constant. Each of these phenomena leads to 
a smearing of the ion-exchange phase transitions trans- 
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forming them into crossovers. The goal of this section is 
to estimate the relative sharpness of these crossovers. 

Consider first the effect of the finite length (still ne- 
glecting the field escape). Close to the first order phase 
transition the free energy admits two competing minima 
(typically at g = and q ~ 1/2) with the free energies 
FQ^ai) and ^1/2(01), see Fig.[7|(we focus on the alternat- 
ing dopants example of section llllj) . Being an extensive 
quantity, the free energy is proportional to the channel 
length: oc L, where h = 0,1/2. Each of these two 
minima is characterized by a certain ion concentration 
nb(ai) = —ai/{kBTL)dFh/dai. In the vicinity of the 
phase transition at a = ac the difference of the two free 
energies may be written as: 



Fo(ai) - i^i/2(ai) 



ai 



OLc 



(33) 



where Anj^j-^ = no(ctc) ~"'1/2(Q!c) is the latent concentra- 
tion of ions across the transition. Taking a weighted sum 
of the two states, one finds that the concentration change 
across the transition is given by the "Fermi function" : 



An(ai) 



Ar 



(34) 



where AA^ = Anj^j^L is the total latent amount of ions 
in the finite length channel. This gives for the transi- 
tion width (ai — ac)/ac oc 1/AA^. Therefore the transi- 
tion is relatively sharp as long as AA^ ^ 1. For small 
enough 7 the number of ions entering or leaving the 
channel at the phase transition is almost equal to the 
number of dopants: AA^ ^ ^dopants ~ iL/xt- The 
necessary condition of having a sharp transition, there- 
fore, is to have many dopants inside the channel. For 



example, for the transition of Fig. ^1 Ar 



0.87/xT 



and 7 = 0.1, one finds AiV « O.OSL/xt- At larger 7 
the number of dopants increases (for fixed length i), 
but the relative latent concentration is rapidly decreas- 
ing, see Fig. ^] Employing Eq. @, one estimates 
AA^ « ^{\ — 2er^/'^^^'>)L / XT this expression is maximized 
at 7 ~ 0.15, where AA^ « Q.IL/xt- We notice, in pass- 
ing, that even for AA^ < 1, being plotted as a function 
of log a2 , the function H37(l still looks as a rather sharp 
crossover. 

We turn now to the discussion of the electric field es- 
cape from the channel, which happens due to the finite 
ratio of the dielectric constants of the channel's interior, 
Ki, and exterior, ^2- Fig. 1251 shows the electrostatic po- 
tential of a unit point charge placed in the middle of the 
channel with ki/k2 = 40, 15]. The potential interpo- 
lates between —e/(Kix) at small distances, x ^ a, and 
2C/oo/e — e/K2X, where Uoo = ^ / {kio^) , at large dis- 
tances a; > ^. Here the length scale ^ ~ a^/^ifK^ « 6.8a 
is the characteristic escape length of the electric field 
displacement from the interior of the channel into the 
surrounding media with the smaller dielectric constant. 
The quantity f/oo = Ul{0)2£^/ L is the excess self-energy 
of bringing a unit charge inside the infinite channel. 
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FIG. 25: Electrostatic potential <l>(a;) in units of e/nia of a 
point charge in an infinitely long channel as a function of the 
dimensionless distance x/a along the channel axis. Here a 
is the radius of the cylinder, and ki/k2 = 40. The full line 
is an exact solution of the Laplace equation^*—. The dotted 
lines are 3d Coulomb potentials: the lower one is —e/nix; the 
upper one is — e/K2X, where e^ao = 2Uoo and Uoo is the 
self-energy of a charge in the infinite channel. The dashed 
line corresponds to Eg 1351 



In between the two limits the potential may be well 
approximated by the following phenomenological expres- 
sion (dashed line in Fig. I25|l : 



1 



(35) 



Our previous considerations correspond to the limit ^ — s- 
00 (save for the last term). The last term in Eq. H35() 
originates from the fact that in the immediate vicinity of 
the charge |a;| ^ a the electric field is not disturbed by 
the presence of the channel walls. As a result the length 
of « 1.1a is excluded from paying the excess self-energy 
price. This leads to (typically slight) renormalization of 
the effective concentrations o; — > Op// . The more detailed 
discussion may be found in Ref. p ; below we neglect the 
last term in Eq. H35|l . 

One can repeat the derivation of section IVII with 
the potential Eq. (|35l) . employing the fact that <I>~^ = 
{2Eo)^^S{x — a;')[^^ — d^]. As a result, one arrives at 
Eq. H21|l with the modified Hamiltonian H = {idgY — 
2q!iC0s6' -I- {xt /2£)'^9'^ . Since the last term violates 
the periodicity, the quasi-momentum q is not conserved. 
However, in the limit xy/C ^ 1 one can develop the 
quasi-classical approximation over this small parameter. 
Transforming the Hamiltonian into the momentum rep- 
resentation, one notices that q{x) is a slow quasi-classical 
variable. As a result, the partition function of the chan- 
nel with the finite ^ may be written as: 



Z = y Vq{x) exp l—Jdx 



1 



— (5,<z(x))^ + — e 



XT 



XT 



(36) 
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where Fq = ksTt^q^ L/xt is the free energy as function 
of g in ^ — > c» hmit (no electric field escape). This ex- 
pression shows that there are no true phase transitions 
even if e^q^ possesses two separate minima. Indeed, due 
to its finite rigidity the q{x) field may form domain walls 
and wander between the two. As a result, the first or- 
der transition is transformed into a crossover. Formally 
Eq. (|36|l defines the "quantum mechanics" with the po- 
tential ^ e"q \ The smearing of the transition is equiv- 
alent to the avoiding crossing intersection due to tun- 
nelling between the two minima of the e^'^'' potential. Us- 
ing this analogy, one finds for the concentration change 
across the smeared transition: 



1 



a/ (ai - Q!c)2 -I- a2^2 



, (37) 



where 5 is the WKB tunnelling exponent: 

1/2 



J = exp <( --^ / dq\Jeq^ - e^,"^ 



(38) 







As a reasonable approximation for e^"' (etc) one may 
use (c.f. Eq. glj) ef\ac) = Uc/{8Ul{0)) cos{4TTq), 
where Uc is the transport barrier at the critical point. 
Substituting this expression in Eq. H38|) one estimates 
S « exp{-^/{2TTXT)VUc/UL{0)}. Using f « 6.8 a and 
Uc « 0.2C/l(0) (cf. Fig.EJ, one obtains 6 « exp{-/s/a}. 
Thus for channels with a < Ib one obtains 5 < 0.4 and 
the crossover Eq. (|57|) is relatively sharp. 
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FIG. 26: The electrostatic potential for cations across the 
channel with 5 dopants, considered in the text. 



the bare parabolic barrier of Fig. ^ we arrived at the wide 
well with the almost flat bottom fFig. I26|l . 

Unlike the L-dependent self-energy barrier considered 
above, the role of the Donnan potential is charge sensi- 
tive. In the negatively doped channel the Donnan po- 
tential well exists only for the cations while anions see 
the additional potential barrier. Thus, the Donnan po- 
tential naturally explains why in negatively doped bio- 
logical channels the cation permeability is substantially 
larger than the anion one. 

The contact potential Ud may be augmented by the 
negative surface charge of the lipid membranei^ or by 
affinity of internal walls to a selected ion, due to ion- 
specific short range interactionsi*^. It seems that biolog- 
ical channels have evolved to compensate large electro- 
static barrier by combined effect of Ud and short range 
potentials. Our theory is helpful if one wants to study 
different components of the barrier or modify a channel. 
In narrow artificial nanopores there is no reason for com- 
pensation of electrostatic barrier. In this case, our theory 
may be verified by titration of wall charges. 



IX. CONTACT (DONNAN) POTENTIAL 

Until now we concentrated on the barrier proportional 
to the channel length L (or escape length f). If a ^ 7 
there is an additional, independent on L, contribution to 
the transport barrier. It is related to the large difference 
in cation concentrations inside and outside the channel. 
Corresponding contact (Donnan) potential Uu is created 
by double layers at each end of the channel consisting of 
one or more uncompensated negative dopants and posi- 
tive screening charge near the channel's mouth. 

For 7 <C 1 one finds ^ U^Lil) and the channel 
resistance remains exponentially large. When 7 grows 
the barrier Ul{^) decreases and becomes smaller than 
Ud = —kBTlii{'y/a), which increases with 7. In this 
case the measured resistance may be even smaller than 
the naive geometrical diffusion resistance of the channel. 

Let us, for example, consider a channel with L = 5 
nm, a — 0.7 nm, xt — 0.35 nm at c = 0.1 M (which 
corresponds to a = 0.035) and 7 = 0.3 (5 dopant charges 
in the channel). The bare barrier Ul{0) = 3.5kBT is 
reduced down to Ul{j) = 0.2kBT. At the same time 
Ud = —2.5kBT. Thus due to 5 wall charges, instead of 



X. CONCLUSION 

In this paper we studied the role of wall charges, which 
we call dopants, played in the charge transport through 
ion channels and nanopores. We considered various dis- 
tributions of dopant charges and salt contents of the so- 
lution and showed that for all of them doping reduces 
the electrostatic self-energy barrier for ion transport. 
This conclusion is in qualitative agreement with general 
statement on the role of transitional binding inside the 
channeliS'. In the simplest case of identical monovalent 
dopants and a monovalent salt solution such a reduc- 
tion is monotonous and smooth function of the salt and 
dopant concentrations. The phenomenon is similar to 
the low-temperature Mott insulator-metal transition in 
doped semiconductors. However, due to the inefficiency 
of screening in one-dimensional geometry we arrived at a 
crossover rather than a transition even in infinite channel. 

A remarkable observation of this paper is that the in- 
terplay of the ion entropy and the electrostatic energy 
may lead to true thermodynamic ion-exchange phase 
transitions. A necessary condition for such a transi- 
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tion to take place is the competition between more than 
one possible ground states. This in turn is possible for 
compensated (e.g. alternating) doping or for mixture of 
cations of various valency. The ion-exchange transitions 
are characterized by latent concentrations of ions. In 
other words, upon crossing a critical bulk concentration 
a certain amount of ions is suddenly absorbed or released 
by the channel. The phase transitions also lead to non- 
monotonic dependencies of the activation barrier as a 
function of the ion and dopant concentrations. For sim- 
plicity we restricted ourselves with the periodic arrange- 
ments of dopants. The existence of the phase transitions 
is a generic feature based only on the possibility of having 
more than one ground state with global charge neutral- 
ity. Thus they exist for arbitrary positioned dopants. In 
reality the phase transitions are smeared into relatively 
sharp crossovers due to finite size effects along with the 
finite electric field escape length, ^. 

We have also demonstrated that the doping can make 
the channels selective to one sign of monovalent salt ions 
or to divalent cations. This helps to understand how bio- 
logical K, Na, Ca channels select cations and how Ca/Na 
channel selects Ca versus Na. The surprising fact is that 
Ca^+ ions, which could be expected to have four times 
larger self-energy barrier, actually exhibit the same bar- 
rier as Na"^. This phenomenon is explained by fraction- 
alization of Ca^+ on two unit-charge mobile solitons. 

We study here only very simple models of a channel 
with charged walls. This is the price for many asymp- 
totically exact results. Our results, of course, can not 
replace powerful numerical methods used for description 
of specific biological channels^^. 

In the future this theory may be used in nano- 
engineering projects such as modification of biological 
channels and design of long artificial nanopores. An- 
other possible nano-cngincering application deals with 
the transport of charged polymers through biological or 
artificial channels. A polymer moves slowly and for ions 
its charges may be considered as static. Therefore, for 
thin and stiff polymers in the channel, the charges on 
polymers can play the same role as doping. As a re- 
sult, all the above discussions are directly applicable to 
the case of long charged polymer slowly moving through 
the channel. Changing the polymer one can change the 
dopants density. 

In a more complicated scenario, the polymer can be 
bulky and occupy substantial part of the channel's cross- 
section. Important example of such situation is translo- 
cation of a single stranded DNA molecule through the 
a- Hemolysin channel^®. In this case, the narrow part of 
the channel, immersed in the lipid membrane (/3-barrel) 
can be approximated as an empty cylinder, while DNA 
may be considered as a coaxial cylinder blocking approxi- 
mately a half of the channel cross-section. The dielectric 
constant of DNA is of the same order as one of lipids. 
Thus, the electric field lines of a charge located in the 
water gap between the two lipid cylinders are squeezed 
much stronger than in the empty channel. This may ex- 



plain strong reduction of the ion current in presence of 
the DNA, which is also different for poly-A and poly-C 
DNA^^a. The latter remarkable observation inspires the 
hope that translocation of DNA may be used as a fast 
method of DNA sequencing. We shall discuss the bulky 
polymer situation in a future publication. 
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APPENDIX A: PHASE TRANSITIONS AT 
LARGE DOPANT CONCENTRATION 

In this appendix we discuss some details of the phase 
transitions at jci ~ 1-7, 7c2 ~ 5.3, jc3 ~ 10.8, etc, visible 
in Fig. ^1 The free energy Fq as function of the order 
parameter q for a few values of 7 in the vicinity of 7ci 
is plotted in Fig. 1271 The minimum, initially at 5 = 1/2 
for 7c < 7 < 7ci, splits into two minima symmetrical 
around 1/2. These two gradually move away from each 
other until they rich q = and g = 1, correspondingly, 
for 7 > 7ci. The continuous variation of the absolute 
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FIG. 27: Free energy in units of 10"^°)7l(0) as a function 
of 5 for 7 = 1.7179330, 1.7179337, 1.7179341 and 1.7179346 
(from top to bottom) with the same 02 ~ 5-10"^. The graphs 
are vertically offset for clarity. The ground state continuously 
changes from g = 1/2 to g = 0. 

minimum (as opposed to the discrete a switch between 
the two fixed minima) suggests the second order phase 
transition scenario, taking place at 7ci. The situation is 
more intricate, however. Namely, there are two (and not 
one) very closely spaced second order transitions, situ- 
ated symmetrically around jci point. At the first tran- 
sition the minima depart from g = 1/2 and start moving 
towards q = and q = 1, while at the second one they 
"stick" to g = 0, 1. Therefore, unless one has a very fine 
resolution (in 7 and/or 0^2), the entire behavior looks as 
a single first order transition. 

The Fq functions of Fig. are very well fitted with 
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the following phenomenological expression: 

Fqil) = a(7ci - 7) cos(27rq) + 6cos(47rg) , (Al) 

where a ^ 6 > 0. For any 7 the ground state corre- 
sponds to the value go which minimizes Fq. Therefore 
go is found from: either cos(27rgo) = (7 — 7ci)a/(4fo), or 
sin(27rgo) = 0. The former equation has solutions only 
in the narrow interval 7ci — Ah /a < 7 < 7ci + Ah /a. In 
this interval of dopant concentrations the minima move 
from g = l/2to(7 = 0, 1. The edges of this interval con- 
stitute two second-order phase transitions located in the 
close proximity to each other. Near the first transition 
I go ^ 1/2 1 — ^7 + 46/a — 7ci, while near the second one: 
|go| ~ \/ 1 ~ Ah I a — 7ci. Therefore the critical exponent 



is the mean-field one: /3 = 1/2. This could be anticipated 
for the system with the long-range interactions. 

It is interesting to notice that the first order transition, 
discussed in the main text, may be also well fitted with 
Eq. IjAip but with negative coefficients a and h. The accu- 
racy of our calculations is not sufficient to establish if the 
subsequent transitions are the first order ones or pairs of 
the very closely spaced second order transitions. As far as 
we can see the sequence of the reentrant transitions con- 
tinues at larger dopant concentrations. Notice, however, 
that the difference of the corresponding free energies (and 
thus associated latent concentrations) are exponentially 
small at large 7. 
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